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Abstract 

In this paper we introduce a 3D phenomenological model for shape 
memory behavior, accounting for: martensite reorientation, asymmet- 
ric response of the material to tension/compression, different kinetics 
between forward and reverse phase transformation. We combine two 
modeling approaches using scalar and tensorial internal variables. In- 
deed, we use volume proportions of different configurations of the crys- 
tal lattice (austenite and two variants of martensite) as scalar internal 
variables and the preferred direction of stress-induced martensite as 
tensorial internal variable. Then, we derive evolution equations by 
a generalization of the principle of virtual powers, including micro- 
forces and micromovements responsible for phase transformation. In 
addition, we prescribe an evolution law for phase proportions ensuring 
different evolution laws during forward and reverse transformation of 
the oriented martensite. 



Key words: shape memory, phase transformation, reorientation, internal 
variables 
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1 Introduction 



In the last years shape memory alloys have been deeply investigated, from 
the point of view of modeling, analysis, and computation. Indeed, these 
materials present many important industrial applications (e.g. aeronautical, 
biomedical, structural, and earthquake engineering) due to their character- 
istic of superelasticity and shape memory effect. 

It is known that the shape memory effect is the consequence of a (re- 
versible) martensitic phase transformation between different configurations 
of the crystal lattice in the alloy: from a high symmetric phase, austenite, 
to a lower symmetric configuration, martensite. Austenite is a solid phase 
(present at high temperatue) which can transform in martensite by means 
of a shearing mechanism. When transformation comes from thermal actions 
(lowering temperature) the result is a multi-direction martensite, in which 
variants compensate each other and there is no resulting macroscopic defor- 
mation. On the contrary, when transformation is obtained by loading, ori- 
ented martensite is formed in the stress direction, exhibiting a macroscopic 
deformation. 

In particular, the research have been developed towards the aim of finding 
a flexible phenomenological model. Some reliable models have been proposed 
to predict the response of such materials. Among the others, we focus on 
two models developed in the framework of phase transitions. 

The first, proposed by Fremond (cf., e.g., [7] and some generalizations 
[1], [5]) describes the behavior of shape memory in terms of (local) volume 
proportions of different configurations of the crystal lattice. More precisely, 
the austenite and two variants of martensite are taken into account. Note 
that the average behavior of different configurations is considered as the 
behavior of the equivalent single variant. The resulting model is able to 
describe phase transformation between different configurations. However, 
the model is obtained assuming that the direction of the transformation 
strain (associated to the detwinned martensite) is known. 

The second model we are considering has been proposed in [TU] and 
then generalized in [2]. In this case one internal variable describes the phe- 
nomenon: the transformation strain tensor. In this case, the internal variable 
describes the direction of martensite orientation. In addition, it leads to a 
simple and robust algorithm, based on plasticity-like return map. Thanks 
to this property the model has been used for implementation within finite 
element codes, allowing the simulation of complex SMA devices. However, 
some secondary effects are not included in this second model, as scalar and 
directional informations are tightly interconnected. 

Thus, one could wonder how to get a deeper description of micro-phenomena. 
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possibly combining the main features of the two different approaches. Ac- 
cordingly, the purpose of this paper is to answer to this idea combining the 
two theories describing secondary effects in the phase transitions as well as 
directional information for the transformation strain. Thus, both scalar and 
tensorial internal variables are introduced, accounting for the phase propor- 
tions (assuming that in each point the phases may coexist with different 
proportions) and for the orientation of the transformation strain associated 
to the detwinned martensite. We consider both proportion and direction as 
internal variables and we write evolution equations for both of them. We re- 
call that an attempt in this direction has been performed in [1]. In particular, 
we prescribe an evolution law to capture asymmetric response of the material 
in tension-compression loading. In [3] a ID model has been introduced to 
describe this kind of phenomenon, using an asymmetric energy depending 
on a tensorial variable. However, this model seems to be hard to be ex- 
tended to higher dimensions. On the contrary, our approach (developing an 
asymmetric evolution theory for phase proportions) can apply to any space 
dimension. 



In this section we detail the derivation of the model. We mainly refer to 
the approach proposed by Fremond to describe the behavior of a thermo- 
mechanical system in terms of state and dissipative variables, as well as 
energy and dissipation functionals (see [7]). The main idea consists in as- 
suming that (microscopic) phase transformations are due to micro-forces and 
micro-movements that have to be included in the global energy balance of 
the system (i.e. generalizing the principle of virtual powers). In particular, 
the equations governing the evolution of internal variables are recovered as 
balance equations (as for the momentum balance). 

2.1 The state and dissipative variables 

As it is known, phase transformations in the alloy are due to the phase 
transitions occurring in the microstructure configuration between austenite 
and twinned or detwinned martensite. In particular, detwinning manifests 
itself mainly through a shear strain so that we introduce a symmetric and 
deviatoric strain which appears in presence of the detwinned martensite. The 
(local) volume proportions of austenite and martensite variants is represented 
by phase parameters 



2 The model 



Xa,Xm,Xs ^ [0,1], Xa + Xm + Xs = ^- 




varf ase 
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More precisely, xa stands for austenite, xm for twinned martensite, and xs 
for detwinned martensite. Furthermore, cZ*'' is the direction of the deviatoric 
strain tensor associated to the detwinned martensite with = {Cs is 

the maximum amount for the detwinned martensite). Indeed, the deviatoric 
strain for detwinned martensite is given by Xsd^"^ ■ Then, is the absolute 
temperature, s{u) the (symmetric) linearized strain tensor {u is the vector 
of small displacements as we restrict ourselves to small deformations). Due 
to the internal constraint (12.11) on the phase proportions (coming by the their 
physical meaning) we can restrict ourselves to consider just two independent 
phase variables {xm,Xs) letting 

Xa = 1 - Xm - Xs, 

where 

0<Xm,Xs<1, Xm + Xs<1- 

Finally, let us use the notation for the elastic component of the strain, so 
that it results 

e = e' + Xsd'\ 
Hence, the corresponding deviatoric strain e is 

e := £ - ^tr (e)! 

I being the identity matrix and tr (■) the trace operator. If a is the Cauchy 
stress tensor, the deviatoric stress tensor S is 

S := cr - ^tr (cr)I = cr - a^J-. 

As far as evolution, this is described by dissipative variables XMti Xst, dl^' , 
and V^. These variables are in particular related to micro-velocities in the 
phase transformation. 

Remark 2.1. Let us comment about the choice of state variables. The main 
idea consists in distinguishing between the norm and the direction of the 
inelastic strain. In this way, we are able to describe the presence of a product 
phase, to which a homogenized strain is associated, and a parent phase, in 
which we find only elastic strain. However, in the parent phase, we can also 
distinguish between the presence of twinned martensite and austenite. Thus, 
we get a more complex and free description of the phenomenon with respect 
to the Souza and Fremond models, which could be useful in some situations. 
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2.2 The energy and dissipation functionals 

We introduce the following free energy functional (depending on state vari- 
ables) as a combination of the energies associated to the single variants (com- 
bined with suitable proportions) and by an interaction energy, accounting 
also for internal constraints 



where 



^(£, XM, X5, 0) = ^el + + ^ch + (2.2) 



'^ei=[^ + ^]{tref + ^\\e-xsd"-f (2.3) 



tr 



2 3 

^id = c,{{e-e^)-9\og 

^ch = (1 - Xm - Xs)hA{e) + XMhM{e) + XshsiO) + hd{e) : d 

= Ik{XM,Xs) + hsi\\d'"\\) + ^^nt{XM,Xs)■ 



Here > is the specific heat, A and /i are the Lame constants; hA,hs, Hm, 
are smooth thermal functions whose regularity will be specified later on (to 
ensure compatibility with thermodynamics). The function is the indicator 
function of the convex set K 

K := {(xAf, Xs) e : < Xm, X5 < 1, Xm + X5 < 1}, 

i.e. it is Ik{xm, Xs) = if {xm, Xs) e K, while Ik{xm, Xs) = +oo otherwise. 
The function /^^ forces = Indeed, it is J^^dld*"!!) = if = 6 

and it is -|-oo otherwise. \E'j„t is a (sufficiently) smooth function accounting 
for interaction energy. As a possible choice for the interaction energy \E'j„t, 
we could simply consider 

^intiXM, Xs) = CmsXmXs + (CamXm + CasXs){1 - Xm - Xs) (2.4) 
+CamsXmXs{^ — Xm — Xs), 

where Cms, Cam, Casi Cams are positive constants. 

Remark 2.2. Note that for {xm,Xs) we have introduced a convex constraint 
forcing (xmiXs) £ K. The constraint on cZ*'' is convex w.r.t. to its norm as 

it is = es. 

Now, let us introduce the pseudo-potential of dissipation, which is a pos- 
itive convex functional depending on dissipative variables, vanishing for van- 
ishing dissipation (cf. [9]). We have 

<f>{XMt,Xst,d';,V9) = \xMt\+Mxs,cT,xst) + Xs\\dT\\ + ^\Ve\^. (2.5) [pseudo 
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Note that, cj) is considered to possibly ensure an asymmetric behavior in 
tension and compression and for forward and backward transformation. This 
is due to the choice of the function 05 (noting that it possibly depends on 
the stress and xs)- Indeed, this choice generalizes the classical situation for 
rate-independent systems, where it is 

05(X5,cr,X5t) = |X5t|- (2.6) 

Actually, (j)s is required to be rate independent with respect to xst- Hence, 
accounting for a possible dependence in the evolution on the stress (e.g., 
for tension-compression behavior) and on the volume of already detwinned 
martensite, we get further possible example 



<^5(X5, cr, Xst) = d{xs, o-){xst) + \xst\ (2.7) specifico 



where (/)"*" = / if / > and (/)^ = if / < and is a sufficiently smooth 
function. From now on we deal in particular with (12. 7p . 

asim Remark 2.3. Note that we could refine the model, e.g. adding in (12. 5p a term 
as d{xsy (^){Xst)~ for decreasing evolution of the product phase. 

2.3 The equations 

We consider a smooth bounded domain Q C with F = dQ split into 
Fi U F2 (with Fj disjoint subset, Fi with strictly positive measure). 

We assume that a generalized version of the principle of virtual powers 
holds, accounting for internal microforces responsible for phase transitions. 
Thus, the first principle of thermodynamics reads as follows 

et + div q = r + a: St + BmXmi + BsXst + B : in 1] (2.8) 

the right hand side being the power of interior forces and the heat source 
r. Here, e is the internal energy, q the heat fiux, {Bm,Bs) and B internal 
(microscopic) forces responsible for the phase transformation (i.e. the evolu- 
tion of internal variables). The heat fiux satisfies boundary condition {h is a 
known fiux through the boundary) 

q-n=/i onF. (2.9) 

Hence, by the principle of virtual powers we get the quasi-static momen- 
tum balance 

-divcr = f in n (2.10) 



enl 
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with boundary condition 



u = on Fi (2-11) boundul 

o-n = t on T2 (2.12) 



boundu2 



f being a volume force, while t is a traction applied on a part of the boundary. 

Analogously, the evolution of the phases depends on internal forces which 
are included in the energy balance of the system. Thus, we get two bal- 
ance equations, one for the evolution of the phase proportions (related to 
{Bm, Bs)) and one for the evolution of the tensor d*^ (related to B), i.e. 

{Bm,Bs) = (0,0) inn (2.13) 
B = inn. (2.14) 

2.4 The constitutive relations 

We need to prescribe constitutive relations for the involved physical quanti- 
ties. The internal energy is 

6 = "^ -6rj 
where the entropy r] is prescribed by 

v = -^ = cshge- h'^ie)ii - XM - xs) (2.15) 
- h'j^{e)xM - h's{e)xs - Kio) 

The Cauchy stress tensor is 

cr = S + (7^1, 

with 

c^tr 6X3 



dm A 2 , 

0"m. = TT = + 77/^ tr £, 



and 

S = 

de 



dm , , 

S = — = 2//(e - xsd )• (2.16) I Sgrande 



Hence, we get 



(5,„ Bs) = {Bll, Bf) + {Bt„ Bi) = — ^ + J\ (2.17) 

0[Xm,Xs> o[XMt,Xst> 

More precisely, letting 

s{x) = a X ^ 0, s(0) = [-l,l], 
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and 



there holds 



H{x) = 1 if x > 0, H{x) = if X < 0, H{0) = [0, 1] 



B'lf = -hAie) + hM{9) + ^ + 7M 



int 



dXM 



and (choosing (ps as in (12.71) ) 



(2.18) scelgoBl 



B 



d _ cxPs 



St 



dxs 

s{xst) + d{xs:0-)H{xst) 



(2.19) scelgoB2 



with 



Finally, we consider 



(7m,7s) e dlKiXM,Xs)- 



(2.20) 



(2.21) 



where 



B'^ = Xss{dr) 



(2.22) 



using the notation 



and letting 



s«) = if ^ 0, s(0) = {w : ||w|| < 1}. 



ieydkX\\d''-\\) = dkX\\d''-\ 



As far as the heat flux, we assume (Fourier law) 



(2.23) 
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3 The PDE system 



3.1 The first principle 

Combining constitutive relations with the balance laws, we get the PDE 
system we deal with. First let us discuss the energy balance, from which 
we show that the model is thermodynamically consistent. The equation 
governing the evolution of the temperature is recovered from fl2.8p . After 
using the chain rule and by the constitutive relations, we get 

from which the second principle of thermodynamics follows, once ^ > 
(it is the absolute temperature). Note in particular that we have strongly 
exploited the fact that 9$ turns out to be a maximal monotone operator 
with G (?$(0). The resulting equation is 

Otics - e{hi,{e)xM + h"s{e)xs + h\{e){i - xm - xs) + h^'{0) ■■ d'l) 

(3.25) 

+ eh'j,{e){xM + xs)t - oh',,{e)xMt - eh's{e)xst - eh^{e) ■. df - m 

= \XMt\ + \xst\ + d{xs^(^){.Xst)''\xst\ +Xs\\(i^[\\- 

In particular, we have to assume that hA,hM,hs,hd are smooth functions 
such that 

(c, - 9{h'i,i9)xM + h"s{9)xs + -XM- Xs) + hj'{e) : d*^)) > C > 0. 

3.2 The evolution 

Combining constitutive relations with momentum balance, it follows 

- div ((A + \ii)ti el + 2/x(e - Xsd'')) = f , (3.26) 

combined with (12.111) . (I2.12p . Then, by definition of Bm and Bs, the evolu- 
tion equations for {xMt,Xst) are written as 



siXMt) + {hM{9) - hA{e)) + ^^P^ + = (3.27) I evlphase 

OXm ' 



and 



s{xst)+d{xs,(T)H{xst) + ihsiS)-hAm-Me-Xsdn : d*^+75 = (3.28) ev2phase 
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where 

hM,ls) e dlxiXAuXs)- 
Note that OIk^xm, Xs) = (0, 0) if {xm, Xs) belongs to the interior of K, while 
it is given by the normal cone to the boundary if {xm, Xs) ^ 9K. 
Finally, the evolution equation for d*^ is given by 

X5s(cZf ) - 2/iX5(e - xsd''') + h,{e) (3.29) 



dxs 



7d*^ = 0, 7ea/e.(||rf 



tr I 



Remark 3.1. Note that the coefficient xs of the evolution term s(d/) en- 
sures that in the absence of detwinned martensite there is not dissipative 
contributions involving d*''. 



ev3phase 



3.3 An equivalent formulation 

Let us now introduce {Bm and Bs are defined as in fl2.18p and f l2.19p ) 



FsiBs 



nd\ 



I Tjrid] 
\Bm\ 

ynd 



Bf\-R{Bf,xs.<T 



(3.30) 
(3.31) 



where 



R{B^\xs.<T) = lii Bf<Q 



and R{Bf, xs. o") = 1 + d{xs. <t) if Bf > 0. 
Then, we can rewrite the evolution of the phases (13.271) and (I3.28P as follows 



XMt — C 



A/- 



TDnd 
Dndl 

I -"m I 
Bs 

I R"'^! ' 



Xst — Cs 

C^F^ = 



F,{B^') < 0, 



M,S, 



M, S. 



(3.32) 

(3.33) 

(3.34) 
(3.35) 



Note that play the role of yield functions (see, e.g., [8]). 
Analogously we may introduce 

tnd\ 



Xs 



letting 



with 



'tr 



d 



CdFd = 0, 



B 



nd 



|gn.(i|| ' 

Fd<0. 
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4 Some examples 



In the following we explore the model performances limiting the discussion 
only to the case of a proportional loading state, i.e., neglecting the reorien- 
tation process. Accordingly, to simplify the discussion, we set h^^O) = 0, 
d*'' in the same direction of S (and e), assuming also = 0. Under these 
simplifying positions, we may set (see 02.161) ) 

B""^ e d*" , , 
(4-36) 



|B"'^|| Hell lld*^ 



and 



||S|| = 2MI|e||-X5). (4.37) 

Moreover, we distinguish between two different possible situations, one 
in which we consider only evolution of the stress-induced martensite and one 
in which we consider only evolution of the temperature-induced martensite, 
as discussed in the following. For both problems we start from a material 
completely in austenite (i.e., xs = Xm = 0). 



4.1 Case 1: temperature- induced effect 



For this problem we assume to start from a = and to vary only the tem- 
perature. Accordingly, only a variation of xm can be produced. 

The problem is governed by the following set of equations (see fl2.18p ): 



TDnd 



-hA{9) + hM{e) 



int 



+ 7m, 



M 



Fm{B 

XMt = 

. CmFm 



nd\ 
M ) 



Cm 



= \Blf\-l 

TDnd 
Bm 

I n^d, \ ' 

nd\ 



(4.38) 



Fm{BZ)<0. 



We assume to first properly cool and then heat the material (see Figure ??). 

In Figures ??-?? we report the evolution of the thermodynamic force B^f 
versus the temperature 9, of the temperature-induced martensite xs versus 
the temperature 9, of the quantity 7m versus the the temperature 9. 

It can be observed that during cooling the model is able to reproduce 
a process in which the multi-variant martensite is produced and then dur- 
ing heating a process in which the multi-variant martensite is progressively 
extinguished. The forward and reverse phase transformations are perfectly 
symmetric. 
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4.2 Case 2: stress-induced effect 



For this problem we assume to start from a = and to vary only the stress. 
Accordingly, only a variation of xs can be produced. 

The problem is governed by the following set of equations ( (12. 191) ): 



= -hA{e) + hs{d) + 



2/^(1 |e| I - Xs) + 7S, 



Fs{B\ 

Xst = Csj 
I CsFs = 



Bf\-R{Br,Xs,o- 



nd 



jDnd 
T}nd \ ' 



(4.39) 



FsiB'f) < 



where (letting (12. 7p holds) 

R{Bf,Xs,cT) 
RiBf,xs,(T) 



1 if Bf < and 



d{xs,(T) if i?g^>0. 



(4.40) 



We assume to first properly load and then unload the material (see Figure 
??). 

In Figures ??-?? we report the evolution of the thermodynamic force B^'^ 
versus the applied stress a, of the stress-induced martensite xs versus the 
applied stress a, of the quantity 75' versus the applied stress a, of the applied 
stress a versus the strain e. 

It can be observed that during loading the model is able to reproduce a 
process in which the single-variant martensite is produced and then during 
unloading a process in which the single-variant martensite is progressively 
extinguished. The forward and reverse phase transformation are unsymmet- 
ric. 
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